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Abstract. We propose a fast and efficient technique to create classes of highly 
entangled states of trapped ions, such as arbitrary Dieke states and superpositions 
of them, e.g. NOON states. The ions are initialized in the phonon ground state and 
are addressed globally with a composite pulse that is resonant with the first motional 
sideband. The technique operates on comparatively short time scales, as resonant 
interactions allow one to use the minimum laser pulse area. The number of single pulses 
from the composite sequence is equal to the number of ions, thus the implementation 
complexity grows only linearly with the size of the system. The approach does not 
require individual addressing of the ions in the trap and can be applied both inside 
and outside the Lamb-Dicke regime. 
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1. Introduction 

Entanglement is the most distinctive feature of quantum states involving many particles. 
Within the framework of quantum information science, it may be viewed as a 
resource for the processing of information in ways not permitted by classical logic 
[1]. Entanglement has various physical applications such as dense coding, quantum 
teleportation, quantum cryptography, quantum metrology, etc. [2], which are essential 
for quantum communication and information processing (QIP). It is indisputable that 
entanglement plays a key role in QIP as quantum computers are implemented by many- 
body systems, generally characterized by multi-partite entangled states. Since the 
primary resource for quantum computation is a Hilbert-space dimension, which grows 
exponentially with the available physical resources [3] , the benefits of a quantum over a 
classical computation increase with the size of the physical system. This has inspired an 
intensive research aiming to create and study the properties of multi-partite entangled 
states. 

A very prominent class of such states are Dicke states originally introduced 

in [4]. They contain a given number of excitations n (qubits in state |1)) shared evenly 
amongst all qubits: 



where Pk denote the set of all permutations of the excitations and C„ = 
N\/ [n!(A^ — n)\]. Notably, Dicke states are immune against collective dephasing, which 
is a dominant source of decoherence in various systems, such as trapped ions [5]. 
Therefore, while still offering exponential dimensionality (of C^^2 ~ 2^/ a/ 
example for n = N/2), the Dicke manifold can be used as a decoherence- free 
computational subspace, as in [6]. Dicke states generalize W states, which can be used 
for quantum communication [7]. Furthermore, Dicke states exhibit genuine multi-partite 
entanglement [8] , which is robust against particle loss and is highly resilient vs external 
perturbations and measurements on single qubits [9]. Thus, Dicke states can serve 
as a versatile resource for the preparation of multiparticle entangled states; through 
projective measurements on some of the qubits one can obtain states from various 
entanglement classes. Due to their robust entanglement, these states are particularly 
well suited for the experimental examination of multi-partite entanglement and can be 
used to test fundamental concepts of quantum mechanics. 

Theoretical proposals exist for the generation of Dicke states in a number of physical 
systems, including ensembles of neutral atoms [10], trapped ions [11, 12], quantum dots 
[13] and using linear optics [14]. Of these ion traps are perhaps the best suited for 
their unparalleled level of experimental control. We notice, however, that the existing 
trapped-ion proposals possess one or a combination of the following drawbacks: (1) they 
cannot create arbitrary but only particular Dicke states, and thus do not offer a general 
approach; (2) individual ion addressing is required, which poses significant experimental 
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challenges to scalability; (3) the number of physical interactions needed scales very fast 
with the system size; (4) initialization in a particular Fock state is required; (5) Dicke 
states are achieved with some probability and post selection is required; (6) adiabatic 
techniques require in general very long interaction times. Consequently, extending the 
proposed techniques for larger system sizes remains a formidable challenge. 

In this paper, we take a different approach and propose a simple, general and very 
efficient technique for the creation of a large class of highly entangled states in systems 
of trapped ions. These can be arbitrary Dicke states and any superposition of these, 
such as NOON states, which are invariant, up to a phase, under the exchange of any 
two ions. Besides its generality, the proposed technique is particularly advantageous 
due to several features. It uses composite pulse sequences - a series of laser pulses, 
each with a particular area and phase. Composite pulses are a conceptually simple 
and very powerful control tool, which enjoys large popularity in experimental physics. 
Their simplicity and efficacy in controlling quantum systems stem from the basic 
physical notion of interference. Thus our technique requires much fewer interaction 
steps compared to the traditional approaches, exploiting quantum circuits of a vast 
number of concatenated one- and two-qubit gates: the number of pulses in our approach 
is equal to the number of ions. Therefore it grows only linearly with the system size, 
thus offering only moderate levels of experimental complexity. Another advantage of 
our method is that it assumes collective interaction with all ions and does not require 
to manipulate exclusively individual ions or pairs of such with focused laser beams, 
which often presents a principal experimental challenge. The laser fields are resonant 
with the first motional sideband transition of the ions, which results in short interaction 
times, as opposed to adiabatic techniques. Our technique is apphcable also outside the 
Lamb-Dicke regime. This offers the potential to overcome various detrimental effects, 
such as light shifts and off- resonant excitations [15], which might occur in experimental 
implement at ions . 

2. Model 

2.1. Hamiltonian 

We consider ions confined in a linear Paul trap, which are cooled to their vibrational 
ground state. Each ion has two relevant internal states |0) and with respective 
transition frequency ujq. The linear ion crystal interacts uniformly with a laser pulse 
tuned on one of the sidebands of the (longitudinal) center-of-mass mode, with frequency 
cul = tJo i i^tr, where cutr is the axial trap frequency. The plus sign stands for the blue 
sideband, while the minus sign is for the red sideband. After making the optical and 
vibrational rotating-wave approximations, the interaction Hamiltonian in the interaction 
representation for the red and blue sidebands, respectively, has the form [16] 




-hgit) ya{r^)J+ + a}{r,)J. 



(2a) 
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Hi,B = ^hgit) [a^iv)J+ + (^iv)J~] , (26) 

where g{t) = ?7f2(t) exp(— ?7^/2)/-\/iV is the couphng of the internal atomic states to 
the vibrational mode, producing pulse area A = J^^g(t)dt, rj is the single-ion Lamb- 
Dicke parameter, Q{t) is the real- valued time-dependent Rabi frequency. Here a (77) = 
Y2'^=o(^~^ ^niv'^)\^) {''^'^ M (^Kv) = [^(^7)]^ ^he phonon lowering and raising 
operators, L'^{x) being the generalized Laguerre polynomial. J"^ = J2k=i where 
(fk is the phase of the laser field interacting with the kth ion, and = |lfc)(Ofc| and 
cr^ = |Ofc)(lfc| are the raising and lowering operators for the internal states of the kth 
ion. In the Lamb-Dicke limit the operators 0^(77) and 0(77) become the ordinary creation 
and annihilation operators of the center-of-mass phonons. 

We now perform the transformation \lk) — )■ \lk)e~^'^'', thereby incorporating the 
laser phase into the atomic states. As a result, the Hamiltonian is recast in terms of the 
usual pseudospin operators = ^^i^^- The energy pattern splits into manifolds 
corresponding to n atomic and u motional excitations. Fig. 1(a). For red-sideband 
interaction we have n + u = ttir, i.e. the total number of quanta is conserved, while 
for blue sideband we have n — v = m-Q, i.e. the difference of the quanta is conserved. To 
create symmetric entangled states a suitable choice is = N and m-Q = 0. Though 
our method is equally applicable for both, in what follows we will assume blue-sideband 
interaction. 

Figure 1(a) depicts all states, which are accessible if one starts from the ground 
state |000)|0). The example given is for a chain of = 3 ions. The ions interact with a 
blue-sideband laser field, uj^ = uq + Utr, which couples equally each ion's internal state 
to the collective motional center-of-mass mode: |Ofc)|z/) ^ |lfc)|z/+ 1). The system is 
described by the anti-Jaynes-Cummings Hamiltonian {2b). To this end, we adopt the 
wavefunction notation lip)]^), where \ip) = \q1q2 ■ ■ ■ Qn+i) is the collective internal state 
of the ion qubits, with g/c = or 1, and |z/) is the vibrational Fock state of u phonons. 



2.2. Hilbert space factorization 

In order to study the dynamics of our system, it is convenient to introduce a new basis, 
which consists of the set of eigenvectors of the two commuting pseudospin operators 
and Jz, where = ^(J+J_ + J~J+) + Each state is assigned two quantum 
numbers, j and m^, respectively. Since commutes with the Hamiltonian, the 
Hilbert space factorizes into a set of decoupled chains with different values of j; the 
Hamiltonian preserves j. The meaning of this becomes more transparent if we notice 
that P = J2k,i=i l^ki - i(-l)^'='l, with Ski = a+af + a^a+ + ia^.^a,,, + ^-l)-^^'! 
denoting the action of the swapping operator, which exchanges ions k and /. Therefore, 
each j stands for a particular symmetry with respect to exchanging ions. Since our class 
of target states is invariant under the action of each 5*^;, it comprises the eigenstates of 
Ski with unit eigenvalue. Thus the eigenvalue of is (1 + N/2)N/2. Hence, the chain 
containing our states is assigned j = N/2, and by analogy with the traditional angular 
momentum operators, the number of states is equal to 2j ' + 1 = + 1. These are all 
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Figure 1. (a) Linkage pattern of the collective states of iV = 3 trapped ions equally 
coupled to their common center-of-mass mode by a uniform laser field. The difference 
between the number of ionic excitations n and the number of vibrational phonons v = n 
is conserved. The laser beam is tuned to the blue-sideband resonance, wl = i^o + '^tr- 
(b) Chains comprising symmetric Dicke states for iV = 3. These constitute 

the whole set of states that can be accessed by a uniform laser driving starting from 
the state |00 • • • 0) |0) . The states are coupled resonantly on the first blue sideband with 
^v-i.u{t) being the coupling strengths. 



Dicke states |W^^)|z/), with n being the number of atomic excitations: n = 0,. . . ,N. 
Hence, all states we are interested in are contained in a single chain, and are coupled by 
the Hamiltonian (26) in the order given (Fig. 1(b)). As long as the interaction with the 
ions is distributed uniformly, the dynamics is enclosed in this chain as it gets decoupled 
from the rest of the Hilbert space. The states differ by energy, which is measured by 
the operator J^. Its eigenvalues mj vary from — j to j and define the number of excited 
ions, n = j + nij = N/2 + mj. 

For the following analysis, it will be necessary to go further and calculate the 
coupling coefficients in the new basis. The coupling between the neighbors \j,mj) and 
\j,mj — 1) follows immediately from the matrix elements of the operators a^(?7), a{ri) 
and J±: 

Xu^iAt) = g{t)Ll^i{v^)VN-u + l. (3) 

In the following, for conciseness we will consider only operation inside the Lamb- 
Dicke regime, which requires 77 ^ 1. 
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3. Implementation 

Our method begins witli tlie initialization of the string of ions in the collective internal 
and vibrational ground state \ip)\i^) = |00 ■ ■ ■ 0)|0) [17]. For ease of notation, the indices 
in \ip) denoting states will be omitted hereafter. 

3.1. Creation of Dicke states 

The A^-ion Dicke states \Wj^)\u) are constructed in the following way. A series of 
A^ pulses is applied globally on all ions, each pulse having a particular area Ak and 
phase (pk- The pulses are resonant with the first blue sideband relative to the center- 
of-mass (COM) mode, i.e. the carrier frequencies are wl = cjq + i^tr- Thereby the state 
|00---0)|0) is coupled to all states |iy^)|n), shown in Fig. 1 (b). Because resonant 
interactions are employed, the dynamics is defined only by the pulse areas Aj. and does 
not depend on the temporal pulse shape. 

To study the effect of this interaction we derive the propagator U{Ak, (p^) describing 
the dynamics of the chain of Dicke states (Fig. 1(b)), subject to a laser pulse of 
area Aj. and phase 0^. This is done by exact diagonalization and exponentiation of 
the Hamiltonian, U{Ak,(f)k) = exp ^— ■ie*'^'= / -ffj^dtj. The total sequence of A^ pulses, 

having area Atot = XlfcLi ^fc, is represented by the propagator 

f/tot = U{An, <f)N)U{AN.u <pN-i) ■ ■ ■ f/(Ai, 0l). (4) 

We fix (pi = 0, which defines our phase reference. Hence, the total propagator f/tot is 
defined by A^ pulse areas and A^ — 1 phases, a total of 2N — 1 variables, which can be 
varied as free parameters. 

Dicke states, or various superpositions of these, can be obtained for specific sets of 
parameters. The parameters are determined numerically through minimization of the 
distance T>{c,t) between a calculated state |c) = ?7tot|00 • • -0)10) and a target state \t), 
defined as 

2iv 

vic,t) = j2\{p\i\c)-\m, (5) 
p=i 

where \p) enumerate the collective states of the ions, cf. Fig. 1. Alternatively, 

one can use the fidelity measure J^{c,t) = |(c|t)|^. 

In Table 1 we provide examples of pulse sequences, which yield Dicke states \W^) 
for different number of ions A^. We choose n = [A^/2J, [xj being the integer part of x. 

A numerical simulation of the creation of the state \W^) is depicted in Fig. 2. The 
composite areas A^ and phases 0^ used are listed in Table 1 for A^ = 6. A unit fidelity 
is obtained when the composite sequence comprising A^ (i.e. 6) pulses is applied, the 
total pulse area being Atot = '^k=i = 2.097r. 
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Dicke states, 




Atot (v4i, 0i; ^2, 02; Atv, 





3 2.53 (0.369, 0; 0.484, 2.39; 1.682, 2.976) 



4 2.25 (0.797, 0; 0.493, 1.715; 0.778, 0.549; 0.187, 0.077) 

5 2.11 (0.795, 0; 0.278, 0.403; 0.480, 0.075; 0.223, 0.309; 

0.333, 0.915) 

6 2.09 (0.536, 0; 0.294, 1.53; 0.335, 0.987; 0.172, 1.066; 

0.319, 0.802; 0.434, 0.166) 

7 2.15 (0.107, 0; 0.584, 1.694, 0.562, 1.566, 0.497, 1.313, 

0.039, 1.956, 0.158, 1.301, 0.206, 1.847) 

8 2.46 (0.539, 0; 0.216, 0.389; 0.459, 0.098; 0.251, 1.560; 

0.464, 0.816; 0.25, 0.388; 0.25, 2.078; 0.03, 1.607) 

9 3.35 (0.51, 0; 0.234, 0.83; 0.9, 0.304; 0.19, 2.025; 0.352, 

0.164; 0.379, 0.556; 0.358, 0.097; 0.199, 0.239; 
0.231, 0.471) 

10 3.89 (0.621, 0; 0.367, 1.147; 0.097, 0.994; 0.616, 1.709; 

0.113, 0.263; 0.203, 0.661; 0.579, 0.328; 0.223, 0.831; 
0.775, 0.909; 0.292, 0.462) 

Table 1. Exemplary areas Ak and phases (/)fe (in units of tt) for composite pulse 
sequences, which produce A^-ion Dicke states jW^jf'). We choose n — [N/2\, where 
[xj is the integer part of x. The composite sequences are described by the propagator 
(4) and comprise N phased pulses, tuned on the first blue sideband. It is noteworthy 
that for increasing values of N the total pulse area Atot is less than {N/2)tt (we have 
checked that this property holds also for N > 10). 



3.2. Creation of NOON states 

We can also create arbitrary superpositions of the states contained in the chain of 
accessible states, shown in Fig. 1(b). Of particular interest is the possibility to generate 
NOON states |N^), which are another very important class of highly nonclassical 
entangled states. They can be defined as an equal-probability superposition of two Dicke 
states whereby the excitation is contained either in the internal state or in the 

motional state of the ions: |N^) = \Wi:i)\0) + \W^)\N) = |11 ■ ■ ■ 1)|0) + |00 ■ ■ ■0)|iV) (for 
simplicity, normalization constants are omitted throughout). The ions are maximally 
correlated as measuring the state of one ion determines the state of all ions. 

Various applications of NOON states of both fundamental and practical interest 
have been suggested, such as entanglement enhanced metrology and sub-wavelength 
lithography [23,24]. Different schemes for the creation of NOON states in trapped- 
ion systems have been proposed theoretically [23,25] and realized experimentally [26]. 
However, they are subject to requirements for individual ion addressing, limited number 
of ions or long interaction times where adiabaticity is employed. 

Unlike the previous proposals, our method allows us, by using global addressing, to 
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Figure 2. Fidelity of creating the Dicke state \W^) vs the number of pulses applied. 
These pulses constitute the composite sequence, given by the propagator (4). The 
pulse areas Ak and phases 0fc are shown in Table 1 for = 6. The number of pulses 
required to create \W^) is equal to iV = 6, the corresponding total pulse area being 
Atot = 2.097r. The state is obtained when the fidelity reaches the unity level (iteration 
6 in the figure). 



create arbitrary NOON states. As the interaction is resonant, the states are created on 
short time scales. We perform the same manipulation as for the creation of Dicke states 
- a sequence of laser pulses is applied, each addressing globally the chain of ions and 
having a particular area Ak and a relative phase (pk- The pulses are resonant with the 
first blue-sideband transition. For particular pulse sequences one obtains the coherent 
superposition |M/q^)|0) + \W^)\N). If we logically interchange |0) and |0) ^ 
this state corresponds to the NOON state |N^). 

In Table 2 we provide examples of pulse sequences, which yield NOON states |N^) 
for different number of ions A^. A numerical simulation for the creation of |N^) is shown 
in Fig. 3. A unit fidelity is obtained with the application of the composite sequence 
comprising A^ (i.e. 6) pulses, with total pulse area Atot = J2k=i^k = I.SStt. 

3.3. Rate of creation of the target states 

As already mentioned, the proposed technique operates on comparatively short time 
scales. An estimate for the duration of the composite pulse sequences needed to create 
our symmetric entangled states is given by T^ot = ^tot/fi'- In order to suppress the 
excitation of the extraneous phonon modes (other than the center-of-mass mode), we 
limit the coupling strength from above to g = Wtr/lO [21]. If we assume a typical trap 
frequency of Utr = 4 MHz, as in [22], and Atot ~ Stt as obtained in the above examples 
for the creation of \W^) and |N^), we obtain Ttot ~ 15/is. For comparison, in [11] the 
Dicke state is created adiabatically in 400 /is. 

Importantly, as can be seen from Tables 1 and 2, for increasing number of ions A^ the 
total pulse area A^ot stays below {N/2)7i for Dicke states and below (A^/3)7r for NOON 
states. Therefore, the duration Ttot increases only linearly with A^ and is asymptotically 
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NOON states, 


N^) 


N Atot {Ai, A2, (f)2; ■ ■ An, ( 


Pn) 



3 1.60 (0.696, 0; 0.640, 1.511; 0.259, 1.962) 



4 1.63 (0.402, 0; 0.291, 0.151; 0.667, 1.819; 0.271, 1.465) 

5 1.88 (0.494, 0; 0.249, 0.652; 0.651, 1.271; 0.313, 0.806; 

0.175,1.175) 

6 1.83 (0.285, 0; 0.234, 0.220; 0.099, 0.700; 0.673, 1.178; 

0.403, 0.665; 0.136, 1.023) 

7 2.06 (0.278, 0; 0.300, 0.266; 0.338, 0.034; 0.541, 1.895; 

0.277, 2.138; 0.137, 0.662; 0.187, 0.070) 

8 2.33 (0.259, 0; 0.923, 0.209; 0.346, 0.408; 0.428, 1.572; 

0.003, 1.705; 0.204, 1.216; 0.003, 2.11; 0.162, 1.543) 

9 2.46 (0.395, 0; 0.146, 2.556; 0.186, 1.336; 0.237, 1.854; 

0.680, 0.740; 0.452, 1.660; 0.169, 0.862; 0.007, 
0.222; 0.186, 1.555) 

10 2.93 (0.476, 0; 0.239, 1.247; 0.289, 1.380; 0.256, 0.305; 

0.228, 2.021; 0.415, 0.220; 0.388, 0.749; 0.059, 1.718; 
0.529, 1.823; 0.047, 0.861) 



Table 2. Exemplary areas Ak and phases (f>k (in units of tt) for composite pulse 
sequences, which produce A^-ion NOON states |N^). The sequences are described by 
the propagator (4) and comprise N phased pulses, tuned on the first blue sideband. It 
is noteworthy that for increasing values of N the total pulse area Atot is below (A^/3)7r 
(we have checked that this property holds also for N > 10). 

limited by (A^/2)T'^ and {N/3)TTr, respectively, ~ 8/is being the duration of a tt pulse. 
4. Conclusions 

We have proposed a simple and efficient technique for the creation of arbitrary collective 
states of trapped ions, which are symmetric under exchange of any two ions. These can 
be Dicke states and superpositions of these, such as NOON states. The method uses 
dedicated composite sequences of phased resonant pulses tuned on the first red or blue 
motional sideband of the center-of-mass mode. The composite sequences comprise 
pulses, being the number of ions, thus the implementation complexity and duration 
grow only linearly. This is in contrast to other proposals, which require a rapidly 
increasing number of elementary gates, demanding exclusive interaction with single ions 
or pairs of ions. As opposed to previous proposals, the ions are addressed globally, 
thus individual ion access is unnecessary. Due to the resonant type of interaction and 
because the required by our method total pulse area is as low as {N/2)'it, the states are 
created on a comparatively short time scale. The method is applicable also outside the 
Lamb-Dicke regime. 
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Figure 3. Fidelity of creating the NOON state |N^) vs the number of pulses applied. 
These pulses constitute the composite sequence, given by the propagator (4) . The pulse 
areas Ak and the phases (pk are shown in Table 2 for = 6. The number of pulses 
required to create |N^) is equal to = 6, with the total pulse area being ^tot = I-SStt. 
The desired state is obtained when the fidelity reaches the unity level (iteration 6 in 
the figure). 
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